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; Future surveys for transiting extrasolar planets, including the space-based 

mission Kepler i U, are expected to detect hundreds of Jovian mass planets 
^ '. and tens of terrestrial mass planets. For many of these newly discovered plan- 

^ ! ets, the intervals between successive transits will be measured with an accu- 

racy of 0.1-100 minutes. We show that these timing measurements will allow 
for the detection of additional planets in the system (not necessarily transit- 
ing), via their gravitational interaction with the transiting planet. The transit 
time variations depend on the mass of the additional planet, and in some cases 
Earth-mass planets will produce a measurable effect. 

The one hundred or so currently known extrasolar planets (2) have revealed themselves 
through three types of phenomena: 1) short-duration brightness anomalies in gravitational mi- 
crolensing events caused by planets near the lens star ©, 2) reflex motions of the central 
star (as revealed by radial velocity variations in the stellar spectrum or radio pulse arrival 
times) (|H|5|0[2), and 3) variations in the apparent stellar brightness caused by a planetary 
transit (passage of the planet in front of the star) {8j^J_0iJJjJ2). These approaches provide 
complementary information. Gravitational microlensing measurements primarily constrain the 
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ratio of planet mass to stellar mass. Microlensing surveys are sensitive to Earth-mass planets 
in principle but reveal little orbital information about the planets discovered. Radial velocity 
measurements allow for estimates of the orbital period, eccentricity, and minimum mass of the 
planet. With present technology, radial velocity surveys can only detect planets with masses 
greater than about 10 Earth masses (orbiting low mass stars) <il3\\14l l5). Current transit obser- 
vations allow for the determination of the orbital period and planetary radius. Transit surveys, 
particularly if space-based, will be sensitive to planets as small as Mercury's for the smallest 
stars observed (1). Barring significant improvements in the precision of radial velocity mea- 
surements, the measurement of the mass and radius (and thus the density) of an Earth-sized 
extrasolar planet would appear to be out of reach. 

Here we point out that variations in the time interval between transits, produced via gravita- 
tional interactions with additional planets, allow for the orbital period and mass of the additional 
planet to be determined from transit observations alone. In some instances, when two or more 
planets transit the same star, the density of the planets can be determined. This opens the 
possibility of obtaining considerably more information about transiting planets than has been 
previously thought. 

The time interval between successive transits of an unperturbed planet is always the same, 
because the orbital period is constant. However, it has long been known that the presence 
of a third body can produce short-term variations, in addition to the more familiar long-term 
variations, of the period of a binary il6\\17l . The same is true for planetary systems. It has been 
recently noted that the interval between successive transits of the extrasolar planet HD 209458b 
would vary by ±3 s if a second planet (of mass 10^"^ Mq, period 80-day, and eccentricity 
e ~ 0.4) existed in that system (7^). 

Over the course of their orbits the transiting planet and a second planet exchange energy 
and angular momentum as a result of their mutual gravitational interaction. This interaction, 
greatest at each planetary conjunction, results in short-term oscillations of the semimajor axes 
and eccentricities of the planets, which in turn alter the interval between successive transits. We 
first illustrate this effect by considering our solar system. Figure 1 shows the variation in the 
transit intervals of our terrestrial planets, recorded by distant observers located in the orbital 
plane of each planet. The gravitational perturbations among the planets in our solar system lead 
to transit interval variations ranging from tens of seconds, for Mercury, to thousands of seconds, 
for Mars. The variations for Earth and Venus show oscillations with the 583-day Earth- Venus 
synodic period. 

Next we investigate the influence of an additional planet on the transit interval of HD 209458b. 
We integrate the heliocentric equations of motion of the hypothetical two-planet system (I7PI) and 
test that the relative energy error of the system is less than 10^^^ to assure the reliability of our 
results. We assume that the two planets are coplanar, have orbits that are perpendicular to the 
sky plane, and have initially aligned orbital apsides. We include the mutual gravitational inter- 
actions of the planets and terms that account for the general relativistic (GR) influence of the 
central star (20) but neglect the terms for the GR influence of the planetary masses, terms for the 
oblateness of the star, as well as terms for the tidal interaction with the star. During the course 
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of each integration, we iteratively solve for the central transit times of the inner planet. 

Figure 2 displays the interval between the times of successive transits of HD 209458b, with 
each panel showing the results for a different set of orbital parameters of a hypothetical exterior 
planet. The large spikes in the transit interval, evident in the top three panels of figure 2, 
occur near the times when the outer planet reaches its periastron. Hence, the orbital periods 
of both planets can be determined directly. We choose the period and eccentricity of the four 
hypothetical planets such that their pericenter distances are roughly the same. As a result, 
the largest transit interval variations have comparable magnitude. For clarity, we separate the 
light-time effect (21) (due to the varying distance between the star and the observer as the star 
moves with respect to the center of mass among the star and planets) from the dynamical effects 
described above. 

Whether the presence of a given companion planet can be detected from short-term transit 
interval variations alone depends on the difference between the minimum and maximum transit 
interval. The black lines in figures 3 and 4 show those differences for Jupiter-mass (10^^ Mq) 
and Earth-mass (3 x 10~^ Mq) perturbing planets, respectively, as a function of their orbital 
period (or semimajor axis) and eccentricity. In these calculations the star has mass M^, = Mq; 
the transiting planet has the mass of Jupiter and an initial eccentricity and period of ei = 0.01 
and Pi = 3.0 d, respectively. For each perturber eccentricity, a range of perturber periods 
are tested, in increments of 0. 1 d, starting with the minimum perturber periods that ensures 
that the orbits of the two planets do not initially cross. We follow the same numerical procedure 
described earlier. During each integration, simulating 10"^ d, we record the times of transit, from 
which the minimum and maximum intervals between transits were determined. As expected, 
for a given perturber eccentricity the period variation decreases as the perturber period and 
semimajor axis increase. Likewise, for a given perturber period, the variation is greater for 
larger perturber eccentricity. As indicated earlier, the transit interval variation is primarily a 
function of the periastron distance of the perturber. Comparison of the two figures confirms that 
the magnitude of the effect is proportional to the perturber mass. We also note that, for short 
time spans, the magnitude of the effect is independent of the mass of the transiting planet (a 
result that follows from the equivalence of inertial and gravitational mass). 

Analytic estimates of the variation in transit intervals can be found by integration of La- 
grange's equations of planetary motion (22). Given a transiting planet with semimajor axis ai 
and period Pi, and a perturbing planet with semimajor axis 02 (assumed larger than ai), period 
P2, and mass M2, we find 



where = ^2(1-62) • ^^'^ dashed lines in figures 3 and 4 show the estimate given by the 
expression above. Equation (1) was derived by assuming the perturber follows a parabolic orbit 
with a periastron distance of 02(1 — 62). It significantly underestimates the actual variation 
in transit period for small period ratios and is best suited for 62 > 0.3. The transit interval 
variations increase with Pi when the period ratio of the two planets is held fixed. Thus, the 
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detection of companions is easier for systems where the transiting planet is farther from the 
star. The timing variations are also larger for planets orbiting less massive stars, for a given 
perturbing planet mass. 

As suggested by figure 2, the eccentricity of the outer planet can be estimated from the 
relative magnitudes of the variations Atmax/Atmm- If the period of the perturbing planet is 
much greater that of the transiting planet, the final factor of equation (1) can be ignored. The 
resulting equation can be rearranged to provide an estimate of the mass of the perturbing planet 




given an estimate of 62. 

Large excursions in transit interval variation in figures 3 and 4 occur near small integer ratios 
of the orbital periods of the two planets. These correspond to mean-motion resonances, near 
which the planets undergo larger oscillations of semimajor axis and eccentricity (22). The width 

1 /2 

in semimajor axis of each resonance region is proportion to 02 (M*/M2) and grows rapidly 
with increasing eccentricity (22). The ranges of perturber periods in figures 3 and 4 in which 
the transit interval variations appear irregular correspond to dynamical chaos resulting from the 
overlap of adjacent mean-motion resonances We have excluded the ranges of perturber 
period in which this chaos results in short-term dynamical instability (since such planetary 
configurations are unlikely to be found). 

The two planets of the GJ 876 system provide an excellent example of a 2:1 mean-motion 
resonance in an extrasolar planetary system. Due to the resonant gravitational interaction be- 
tween the two planets, the orbital periods of the inner and outer planet vary from 30. 1 d to 
3 1 . 1 d and 60.0 d to 59. 1 d, respectively, over a libration period of 550 d (7^. Although careful 
photometric monitoring has excluded the possibility of transits of the inner planet of the GJ 876 
system (25), such transit interval variations would be clearly seen if the inner planet transited. 

In some cases, the libration period of a resonant system can be too long for the transit 
interval variations due to the resonance to be easily observed, even if the amplitude of the 
variation is large. GJ 876's short libration period, Pm,, results from the combination of its low- 
mass star (M^ = 0.4 Mq) and massive planets (Mi = 2 Mj^p, M2 = 4 Mjup), since Put, oc 
Pi {MjM^f^ mm (and depends upon the order of the resonance). Note that the 10 d 
integrations used for figures (3) and (4) are long enough to sample a full libration period for the 
low-order resonances. A corresponding system with a solar- mass star, a 1 -year orbital period for 
the inner planet, and jupiter-mass planets in the 2: 1 resonance would have a libration period of 
roughly 50 years. For earth-mass planets the libration period would be nearly 1000 years. These 
long-term effects would not be observable. However, the smaller transit interval variations that 
occur on the time scale of the orbital periods of the two planets, that are superimposed on longer 
period variations, could be observed. 

The feasibility of this technique to detect additional planets in transiting systems depends 
on the physical and orbital properties of the multiple planet systems, as well as the accuracy 
with which the times of transit can be measured. Our expectations for systems with multiple 
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Neptune-mass to Jupiter- mass planets are guided by recent discoveries, as well as by the giant 
planets in our solar system. Fourteen (of 1 17) systems are known to have two or more planets. 
These have period ratios ranging from 2 to 150, and three are known or thought to be in mean 
motion resonances. In the solar system, neighboring giant planets have period ratios of two 
to three. The eccentricities of the planets in the solar system are small, but in many of the 
extrasolar planetary systems the eccentricities are substantial, ranging up to e = 0.7. 

Both the theory of terrestrial planet formation and available observations suggest that the 
typical terrestrial planet system will have a configuration that results in variations in the transit 
interval of hundreds of seconds. Because the escape velocity from the surface of a terrestrial 
planet is smaller than the escape velocity from a solar mass star at roughly one AU, all the solid 
material in the region from a few tenths to one or two AU must either accrete into planets or fall 
onto the star. The resulting planets are as closely spaced as dynamical stability permits (126112 7I) . 
As observed in our solar system and in the pulsar planet system, the period ratios are of order 
two. 

Assuming the observations during the ingress and egress of a transit are well-sampled, the 
error Ot^ (due to photon statistics alone) associated with the measured time of the center of a 
transit of duration is given by 

^-^ ^ {Vtr V p-''\ (3) 
It 

where V is the photon count rate of the star and p = Rp/ is the ratio of the planet radius 
to the stellar radius. Kepler (i), a NASA Discovery mission, will monitor 100,000 A-M dwarf 
stars brighter than apparent magnitude V = 14, looking for transits of detecting Earth-sized 
planets in the "habitable zone" {28). For Kepler, with its 0.95 m diameter aperture, F = 7.8 x 
10^ ]^Q-o.4(y-i4) fQj- a Qf apparent magnitude V. For a Jupiter-sized planet in a 1-year 
orbit about a solar-mass star with V = 14, we find at^ ~ 20 s (p ~ 0.1, t-r ~ 13hr). For 
an Earth-sized planet, at^ ~ 500 s (p ~ 0.01). These accuracies suggest that the transit period 
variations due to the gravitational influence of Earth-mass planets with small period ratios can 
be detected by Kepler. For brighter stars observed with large-aperture ground-based telescopes, 
the presence of more distant additional planets can be detected. For example, observing aV = 9 
star with a 6.5 m aperture telescope, at^ ~ 0.2 s. 

One would like to measure the density of a transiting planet to determine if it is a rocky 
terrestrial planet such as the Earth, an ice/water giant such as Neptune, or a gas giant planet 
such as Jupiter. However, when a transiting planet is perturbed by another planet that does not 
transit, photometry yields an estimate of the radius of the transiting planet but not its mass, and 
the transit timings yield an estimate of the mass of the perturbing planet but not its radius. In 
the cases in which both planets transit their star the masses and radii of both planets can be 
estimated, allowing density determinations. This means that for some subset of the discovered 
planets (double transit systems), radial velocity observations will not be necessary to determine 
their masses and densities. This has important implications for the search for habitable planets. 

If transits of one planet are seen, what is the probability of seeing transits of a second planet? 
Assuming such a planet exists and that the orbits of the two planets are coplanar, the probability 
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that the second planet also transits is simply 01/02 if 02 > oi (if 02 < 0,1 transits of the second 
planet are assured in the case of coplanar orbits). If the orbits of the two planets are mutually 
inclined and one planet transits the center of the star, the probability that the second planet also 
transits is 



where i' is the mutual inclination between the orbits of the two planets rather than the sky plane 
inclination. This assumes sini' > R^/a2, otherwise transits of the second planet are certain. 
For a solar-radius star, a mutual inclination of a few degrees, and a semimajor axis 02 ~ 1 AU, 
the probability is roughly 10%. Thus, such double transiting systems are likely to be found. 

The Kepler team plans to search for transits with a consistent period, depth, and duration (|7]). 
We caution that any detection algorithm that looks for evidence of periodic transits in the pho- 
tometry must allow for variations in the transit period due to the perturbations from unseen 
planets. An overly restrictive test for periodicity might reject some of the most interesting and 
informative planetary systems. The GJ876 system discussed above, an extreme case, could 
easily be rejected as only quasiperiodic. 

We also caution that experiments looking for gradual transit period variations due to slow 
secular trends must allow for the possibility of significant short-term variations in the transit 
period. Transit period variations due to orbital precession induced by the general relativistic 
effects of the central star, the oblateness of the star, or the presence of a distant additional 
planet (EPI) or due to the decay of the transiting planet's orbital semimajor axis as a result of 
dissipative tidal interaction between the star and planet SW^ . require many years to elapse before 
they can be detected; without careful monitoring in the interim to determine the short-term 
variations, the results of such experiments could be easily misinterpreted. 
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Fig. 1. Transit times of the terrestrial planets: The variation of the interval between successive 
transits of terrestrial planets, induced by the other planets in the solar system. To guide the eye, 
the solid line connects the times of each transit. 

Fig. 2. Transit times of HDf 209458b: The interval between successive transit centers of 
HD 209458b as a function of time, with each panel showing the results for a different set of 
orbital parameters for the hypothetical second planet. 

Fig. 3. Variations induced by a Jupiter-mass planet: Variations in the interval between succes- 
sive transits of a planet with Pi = 3 d, ei = 0.01, Mi — 10~^Mq, induced by a second planet 
with mass = lO'^M©. 

Fig. 4. Variations induced by an Earth-mass planet: Variations in the interval between succes- 
sive transits of a planet with Pi = 3 d, ei = 0.01, Mi = 10"^Mq, induced by a second planet 
with mass Ma = 3 x 10~^Mq. 
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Transits of Mars 
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Figure 1 : Variation of transit periods of the terrestrial planets, as would be recorded by distant 
observers located in the orbital plane of each planet. 
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HD209458b (Pj = 3.5248 d, = 0.025) 
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Figure 2: Variations of the transit period of HD 209458b, induced by a hypothetical second 
planet. 
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Figure 3: Variations in the interval between successive transits of a planet with Pi = 3 d, 
ei = 0.01, Ml = IQ-'^Mq induced by a second planet with mass M2 = IO^^Mq. 
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Figure 4: Variations in the interval between successive transits of a planet with Pi = 3 d, 
ei = 0.01, Ml = IQ-^Mq induced by a second planet with mass M2 = 3 x lO^^M©. 
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